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Abstract. 

■ Observations of quasars (QSOs) shining through or close to galaxies offer a way to probe the properties of the foreground 
matter through dust extinction and gravitational lensing. In this paper the feasibility of measuring the dust extinction properties 

, is investigated using the backlitting of QSOs. We test our method to search for QSOs affected by intervening extinction, by 

^ ' matching the coordinates in the SDSS QSO DR3 catalogue with the New York University Value-Added Galaxy Catalog. In total, 

164 QSO-galaxy pairs were found with a distance of less than 30 kpc between the galaxy centre and the QSO line-of-sight at 
the galaxy redshift. Investigating the QSO colours with multiband SDSS photometry, two pairs with galaxy redshifts z < 0.08 
were found to be particularly interesting in that the QSOs show evidence of heavy Galactic type extinction with Ry ~ 3.1 at 
very large optical radii in the foreground spiral galaxies. With the available data, it remains inconclusive whether the two pairs 
can be explained as statistical colour outliers, by host extinction or if they provide evidence of dust in the outskirts of spiral 
galaxies. Deeper galaxy catalogues and/or higher resolution follow-up QSO spectra would help in resolving this problem. We 
also analyse five QSOs reported in the literature with spectroscopic absorption features originating from an intervening system. 

■ These systems are at higher redshifts than the other two and we find in most cases significantly lower best fit values of R v . 
The wide range of preferred values of R v found, although affected by substantial uncertainties, already indicates that the dust 

q . properties in other galaxies may be different from the Milky Way. Furthermore, the available data suggests a possible evolution 

'_h ' in the dust properties with redshift, with lower R v at high z- 
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' 1. Introduction through strong lensing is only efficient for very massive galax- 

■ tes, i.e. less suitable for normal spiral galaxies, especially at 
Understanding the size and composition of the dust grains low re dshifts where the probability for strong lensing is low. 
in galaxies provides important astrophysical information. Furthermore, the time variability of QSOs complicates the in- 
Moreover, it is of critical importance to correct observations ter p r etation of colours from different lensed images, as these 
for extinction losses due to dust, e.g. when measuring cosmo- corresp0 nd to different emission times at the source. In this 
logical distances with supernovae. It is often assumed that the paper> a memo d is proposed where the extinction in galaxies 
reddening properties by dust in all galaxies follow the average is investigated using observed colours of single background 
properties of the Milky Way. The purpose of this work is to QSO images shining through or c i ose to foreground galaxies 
examine ways of improving our understanding of the range of (as shown in Fig [JJ ? taking advantage of the homogeneous 
extinction properties. Earlier analyses of dust in distant galax- spe ctral properties of QSOs. With the current limited data we 
ies have made use of e.g. multiply imaged QSOs (e.g. Nadeau cannot dist inguish between host extinction and intervening ex- 
et al. [19 9JJ Fal co et al. ©99 ) or galaxy backlitting (e.g. Keel t inction. However, the method will be useful with future larger 
& White I200TV The method of using multiply imaged QSOs samples of QSO-galaxy pairs 

has been criticised by McGough et al. <|2005|l since it is im- Even though QSOs show small dispersion in their spectral 
plicitly assumed that one of the images is barely affected by proper ties, a population of redder QSOs have been found in 
the dust or that the reddening laws along the two lines-of-sight studies of QSO colours (see e g Richards et aL I200H . These 
are very similar. In addition, a method based on sources found could be explained invoking extinction by dust, either at the 
central engine, in the host galaxy or in a foreground galaxy. 



Send offprint requests to: L. Ostman There are some indications that the reddening is dominated by 
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Fig. 1. One of the QSO-galaxy pairs where the QSO (SDSS 
Jl 14719.89+522923.1), which is indicated by the arrow, is 
seen through the galaxy. The impact parameter is 5 kpc. The 
redshift of the galaxy and the QSO is zc = 0.0475 and zq = 
1.990. 

extinction from the host galaxy and only rarely from an inter- 



vening galaxy. Hopkins et al. ( 2004 1 found that reddened QSOs 
more often have narrow absorption lines originating from the 
QSO redshift than QSOs unaffected by reddening. No correla- 
tion was found between the reddening of QSOs and the pres- 
ence of narrow absorption lines originating from smaller red- 
shifts. In this work, the presence of QSOs affected by interven- 
ing extinction is investigated further by studying pairs of QSOs 
and foreground galaxies at small projected distances, compa- 
rable to galaxy sizes. As such systems would also fall into the 
category of potential gravitational lensing candidates, we also 
examine their expected (achromatic) magnification. 

In Sec. 13 different parametrisations of dust attenuation are 
discussed. In Sec. [3] our method for finding QSOs affected by 
dust extinction is presented and in Sec. 0] the reddening of a 
sample of QSOs with foreground extinction features are inves- 
tigated. A discussion follows in Sec.[5]and a summary of the 
results in Sec.|6| 

Throughout the paper, the cosmological parameters have 
been assumed to be Qm - 0.3, £2a = 0.7, and Hq - 
70km/s/Mpc. 

2. Reddening by dust 

Within the commonly used extinction laws, the specific dust 
properties can be quantified by the dimensionless reddening 
parameter, the total-to-selective extinction ratio Ry, defined by 



Ry = 



E(B-V) 



(1) 



where Ay is the V-band extinction and E(B - V) is the colour 

excess, 



E(B -V) = (B - V) obs - (B - V) h 



(2) 



(nearly) grey extinction with Ry » 3. Rayleigh scattering 
{A^ °c /L 4 ) would lead to a very steep reddening, Ry ~ 1 .2. 

The typical value found for the reddening parameter in the 
Milky Way is R v ~ 3.1, but it has been suggested that it varies 
between sight-lines in the interval 2 < R v < 6 (see e.g. Draine 
2003 and references therein). In other distant galaxies, values 
of Ry have been derived in the interval 1.5 < Ry < 7.2 (Falco 
et al. 1999 1, where the lower values occur in spiral galaxies and 
the higher values in ellipticals. 

Observations of colours of low -redshift Type la supernovae 
have also been used to constrain the dust properties of their 
host galaxies (e.g. Joeveer 1983 Branch & Tammann 119921 
Riess, Press & Kirshner [T99rj| Reindl et al. 120051 Guy et 
al. [5005 1. Intriguingly, the average reddening parameter has 
been found to be significantly smaller than in the Milky Way. 
However, the effect is hard to disentangle from a possible 
colour-brightness correlation in supernovae or reddening by 
circumstellar dust, possibly significantly different from dust in 
the interstellar medium. Thus, an independent investigation us- 
ing QSO colours could help to break this degeneracy. 

Two parametrisations are often used to describe the extinc- 
tion law in our galaxy: Cardelli, Clayton and Mathis ([H 



The observed colour is denoted by (B - V) b s and the intrin- 
sic colour by (B - V)i ntI . Large grains of dust would produce 



CCM) and by Fitzpatrick ( [19991 FTZ). The difference between 
these parametrisations are most important for large values of 

Ry 

For other galaxies, the extinction curves have only been 
thoroughly investigated for the Magellanic Clouds, where indi- 
vidual stars can be resolved. The Magellanic extinction curves 
are very similar to the Galactic in the visible and the near- 
infrared. In the UV, however, the extinction curve for the Small 
Magellanic Cloud (SMC) differs significantly: the curve is al- 
most linear mA~\ i.e. the Galactic type dust bump at 2175 A 
is not present. The SMC have lower dust and metal abundances 
than the Milky Way and may be more representative of high 
redshift galaxies. There is also evidence that the 2175 A bump 
is not present in damped Lya systems (Pei et al. H999i . Yet, it 
is often assumed that the Galactic laws can be applied to dust 
in other galaxies. 

3. Looking for QSOs affected by intervening 
extinction 

In order to test the feasibility of finding and using QSOs af- 
fected by dust absorption in a foreground galaxy, the coor- 
dinates of the 50748 QSOs in the third data release of the 
Sloan Digital Sky Survey (Schneider et al. 2005 1 were cross- 
correlated with the positions of the galaxies in the New York 
University Value-Added Galaxy Catalog (Blanton et al. 120051 
Abazajian et al. 2004 1. All objects in the galaxy catalogue 
where the classification has been disputable have been re- 
moved. Among these were objects within 0.5 arcseconds from 
a QSO, which often seemed to be multiple detections of the 
same source. We also require that the redshift of the galaxy is 
determined spectroscopically. Only QSOs in the redshift range 
0.5 < z < 3.0 were considered, as many of the QSOs outside 
this range were found to be less accurately described by the 
QSO spectral template made of the HST radio-quiet composite 
spectrum (Telfer et al. 2002 1 combined with the SDSS median 
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composite spectrum (Vanden Berk et al. 2001 ). At higher red- 
shifts, the flux decrement due to the Ly-a break makes the bluer 
bands inefficient. 

In the following, a galaxy will be considered as a fore- 
ground galaxy if the distance from the centre of the galaxy to 
a QSO line-of-sight (the impact parameter) is less than 30 kpc 
at the galaxy redshift. The sample then contains 164 pairs of 
QSOs and galaxies. The smallest impact parameter found with 
respect to the foreground galaxy was 2 kpc. Note that all dis- 
tances are in units of (70km/s/Mpc)///o. Note also that the 
New York University Value-Added Galaxy Catalog is very in- 
complete at z > 0.3. 

For the 164 selected QSOs, the rest- frame B — V colour was 
estimated by K-correcting from the observed bands (corrected 
for Galactic absorption using the values from SDSS given by 
the maps of Schlegel et al. ( 1998 ») using the template QSO 
spectrum. The QSOs with foreground galaxies are only slightly 
redder than the QSOs not associated with a foreground galaxy 
(the mean value of B — V differs by ~ 0.01 mag). The two 
distributions of rest-frame colours were compared using the 
Kolmogorov-Smirnov test, yielding a 78 % probability that the 
difference is only due to statistical fluctuations. Thus, we con- 
clude that most of the QSOs of the selected sample have colours 
that are indistinguishable from the QSOs without associated 
foreground galaxies. Note, however, that this does not exclude 
the possibility that both samples contain QSOs reddened by 
dust in the host galaxy or in intervening galaxies that may or 
may not be resolved. Next, we looked for a possible subset of 
QSOs in the selected sample whose colours were best matched 
with dust extinction at the redshift of the foreground galaxy. 
For that purpose we compared the observed QSO colours with 
synthetic colours derived from the QSO spectral template. The 
template spectrum was redshifted and reddened as if attenuated 
by dust at the redshift of the intervening galaxy. A grid with a 
wide range of values of Ry and E(B - V) was investigated using 
the parametrisation by Fitzpatrick. The expected colours at the 
QSO redshift were used to compute the g, r, i, and z magnitudes 
from the measured w-band magnitude. To find the values of Ry 
and E(B - V) that give the best fit for a certain QSO-galaxy 
pair, we minimise 

X 2 [Rv, E(B - V)] = (X sim - Xobsf^Xsim - Xobs), (3) 

where X slm and X a \, s are the simulated and the observed magni- 
tudes, respectively. *V is the covariance matrix calculated from 
a set of QSOs with similar redshifts (Az < 0.05) without any 
visible foreground galaxies. To avoid a too large impact from 
colour outliers, QSOs with colours differing by more than two 
standard deviations from the mean were removed in the calcu- 
lation of the covariance matrix. 

In the process of selecting which of the 164 QSOs that are 
most likely to be affected by intervening absorption, the fol- 
lowing cuts were made in order to minimise the number of 
"fake" detections while still retaining "good" candidates with 
foreground absorption: 

1 . The fit of the QSO spectrum with the template without dust 
extinction had to be worse (higher x 2 ) than the fit for at 
least 90% of 10000 randomly chosen QSOs without visi- 
ble foreground galaxies. 



2. The best fit for dust extinction must be better than the fit 
for at least 90 % of the randomly chosen QSOs affected by 
a mock foreground galaxy. 

Two QSO-galaxy pairs survived both cuts. The foreground ob- 
ject was a spiral galaxy in both of these systems. For the ran- 
domly chosen QSOs without visible intervening galaxies, only 
0.2 % survived both cuts. Thus, for the QSO-galaxy sample 0.3 
"fake" detections are expected. Using Poisson statistics we thus 
conclude that the likelihood for both events to be false is less 
than 4 % (P(n) = 0.74, 0.22, 0.03, for n = 0, 1,2 respectively). 
However, we cannot exclude the possibility that the two se- 
lected quasars are colour outliers in the poorly understood tail 
of the colour distribution. 

The low number of pairs surviving both cuts consolidates 
our earlier conclusion that the bulk of the sample with 164 
QSOs with foreground galaxies are not significantly different 
from the whole sample. This is also confirmed by a visual 
inspection of all QSO-galaxy pairs showing that, due to the 
relatively large impact parameters, only two of the QSOs are 
clearly behind any visible part of the foreground galaxy. These 
two, however, did not make it through the cuts. Both of the 
QSOs that survived the cuts have an impact parameter close to, 
but clearly outside the visible parts of the foreground galaxy see 
Fig. [2] In both pairs, the surface brightness of the galaxy at the 
position of the QSO is of the order 0.1 % of the central galaxy 
surface brightness. This fact puts limitations to the possible in- 
tervening extinction, unless there are significant amounts of 
dust outside the optical disc of spiral galaxies. In the follow- 
ing, these two pairs will be used to demonstrate our method, 
even though it is possible that the derived extinction is not only 
due to foreground extinction, but could also be due to host ex- 
tinction. 

Information about the two remaining QSO-galaxy pairs, as 
well as the fitted dust parameters are listed in Table ^ The 
best fit values of Ry are listed both for the parametrisation 
by Fitzpatrick and the CCM parametrisation together with the 
error interval under the assumption that the extinction law is 
valid for all positive Ry values. The^f 2 is generally somewhat 
better for the Fitzpatrick parametrisation. In both cases the al- 
lowed values for Ry are consistent with the Galactic average 
(Ry ~ 3.1), as shown in Fig[3]for the CCM parametrisation. 
The parametrisation by Fitzpatrick gives smaller contours (not 
shown). We have also modeled the quasar colours with extinc- 
tion by SMC type dust but the fit is worse than for the CCM 
and FTZ parametrisations. 

As noted in Hopkins et al. (J2>04l, the reddening in the se- 
lected systems could originate in the host galaxy, but such con- 
figurations yield worse fits in both cases. Of course, the redden- 
ing could also be a combination of both foreground and host 
extinction. 

4. Reddening in QSOs shining through galaxies 

Since the New York University Value-Added Galaxy Catalog 
is very incomplete at z > 0.3, we have searched the literature 
for QSO-galaxy pairs with higher galaxy redshifts. Wang et al. 
(2004) have found by a visual inspection three QSO spectra 
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Table 1. Two colour selected QSO-galaxy pairs. The first two columns contain data about the QSO and the following four about 
the galaxy. The positions are given in J2000 and Rp is the petrosian radius in the R band from SDSS. The impact parameters are 
then given both as the angle and the distance at the galaxy redshift. The last two columns contain the best fit values of R v and the 
one cr uncertainty for two different dust parametrisations, Cardelli, Clayton and Mathis ( 1989 ccm) and Fitzpatrick ( 1999] ftz). 
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Fig. 2. The two QSO-galaxy pairs that survived the cuts, SDSS 
J131058.13+010822.2,and SDSS J084957.97+5 10829.0 from 
the top. An arrow points out the QSO in the case where it is not 
clear which object it is. 

in the Sloan Digital Sky Survey (Schneider et al. 12005 1 with a 
2175 A absorption feature at a lower redshift than for the QSO. 
These QSOs are therefore most likely affected by dust in an 
intervening galaxy. The redshift of the intervening system was 
determined using multiple narrow absorption lines. In addition, 
Pindor et al. ( 2003 1 and Johnston et al. ( 2003 1 have found dou- 
bly imaged QSOs where the redshift of the lensing galaxy has 
been determined from absorption features. 

To explore the possibility of dust absorption in these in- 
dividual QSO images, the observed QSO colours were com- 
pared with the simulated colours for different values of Ry 
and E(B - V), in the same way as described in Sec. [3] Note 
that we have not used the differential method that is tradition- 
ally used for determining dust extinction using multiply images 
mentioned in Sec. [2 The best fit values for Ry allowing only 
positive Ry are given in Table|2]for the two different parametri- 
sations. The two last entries in the table correspond to fits us- 
ing two different images of the same QSO. In the following, we 
will use the mean of these two values of Ry. None of the pairs 
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Fig. 3. Confidence levels for the two low-redshift galaxies cor- 
responding to 1 cr (black line), 68 % (dark blue region) and 
90 % (pale blue) as defined by the x 2 test. 

in the table would have survived the cuts in Sec. It is in- 
triguing that the best-fit values of Ry appear to be significantly 
lower than for Galactic dust. 

5. Discussion 

Our QSO colour analysis makes only use of the SDSS QSO 
photometry. However, a visual comparison was also made be- 
tween the candidate spectra and the redshifted template spec- 
trum with and without reddening following the best fit values 
of E(B - V) and Ry. We find that the two candidate spectra 
are consistent with the fit based on the photometric data. These 
comparisons are shown in Fig. |4] However, the spectroscopic 
comparison could not rule out host extinction as opposed to 
extinction in a foreground galaxy. 
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Table 2. Information about QSO-galaxy pairs from the literature. The first two columns contain data about the QSO. The third 
column contains the galaxy redshift and the following four columns contain our best fit values of Ry and E(B - V) for two 
different dust parametrisations and the one <j uncertainty. The last column contains the reference from where the pair was found. 
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Fig. 4. The real spectrum, the redshifted template spectrum and 
the template spectrum after being redshifted and reddened due 
to the best fit dust with the Fitzpatrick parametrisation for the 
two QSO-galaxy pairs. The vertical lines show the average 
wavelengths of the broadband filters. 



We have also looked for narrow absorption lines in the QSO 
spectra. The redshift of the line would tell whether the ab- 
sorption occurs in the host galaxy or in an intervening galaxy. 
Munari and Zwitter ( 1997 1 have found a relation between the 
equivalent width of Na I (5890.0 A) and K I (7699.0 A) and 
E(B - V). This could provide a second test for E(B - V) if 
the lines are detected. Unfortunately, the signal-to-noise of the 
SDSS spectra is not good enough to put limits on E(B - V), 
as shown in Fig. [5] Thus, this is not testable with the public 
SDSS spectroscopic data, but would require follow-up observa- 
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Fig. 5. The SDSS spectrum of the two QSOs with a low- 
redshift galaxy (solid lines) and the redshifted template spec- 
trum with a dip corresponding to the absorption from Na I in 
the intervening galaxy if E(B - V) is given by the best fit value. 



tions with higher spectral resolution and better S/N-ratio. The 
detection of such absorption lines in QSO-galaxy pairs would 
provide evidence for the existence of dust outside the optical 
disc of spiral galaxies, indicative of strong expulsion from dust 
producing regions (Ferrara et al. 1990; Aguirre 1999 1. Without 
such confirmation, the fact that we only have two candidates for 
foreground galaxy extinction that both correspond to positions 
at large optical radii implies that we need to be careful when in- 
terpreting the data. Since the galaxy catalogue is not as deep as 
the QSO catalogue (see Fig.|6j, we expect to suffer from both 
incompleteness in our sample of QSO-galaxy pairs and con- 
tamination in our reference sample. With a deeper galaxy sur- 
vey more QSO-galaxy pairs would be found, also with smaller 
impact parameters. Also, a deeper survey would reduce the 
uncertainty in the calculation of the covariance matrix due 
unresolved foreground objects in the reference QSO sample. 
Forthcoming wide field surveys, e.g. the Large Synoptic Survey 
Telescope (LSST) 1 , could potentially provide over a thousand 
QSO-galaxy pairs where the dust properties in the foreground 
galaxy could be successfully measured as described in Sec.0] 
Systems consisting of a high redshift QSO with a fore- 
ground galaxy close to the line-of-sight are not only interest- 
ing from a dust absorption point of view; it is also a classical 

1 |http : // www . lsst . org 
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Fig. 6. The redshift distribution of the galaxies in the NYU 
VAGC (solid lines) and of the QSOs in DR3 of the SDSS with 
0.5 < z < 3.0 (dotted lines). 

set up for a gravitational lens system. The light from the back- 
ground QSO will become deflected when passing through the 
gravitational field of the foreground galaxy causing images to 
be magnified, distorted and - in rare cases - multiply imaged. 
In order to constrain the magnitude of the lensing effect, we 
model the mass distribution of the foreground galaxies as sin- 
gular isothermal spheres (SIS) with density profiles 



PsisW = 



(T 

2nr 2 



(4) 



The model is characterised by the line-of-sight velocity disper- 
sion <x, which can be estimated from the galaxy luminosity via 
the Faber-Jackson relation (elliptical and lenticular galaxies) or 
the Tully-Fisher relation (spiral galaxies), see e.g. Gunnarsson 
et al. <2005> and references therein. 

For each foreground galaxy, the absolute B-band mag- 
nitude, Mb, is estimated using the publicly available code 
kcorrect v3_2 (Blanton et al. 2003 1. The galaxy type is de- 
termined from the observed galaxy colours and spectra in com- 
bination with visual inspections of the morphologies. 

For the sample of 164 QSO-galaxy pairs with impact pa- 
rameter less than 30 kpc at the galaxy redshift, magnifications 
are generally small (/i < 1.1). The only case where our calcu- 
lations gives a magnification larger than 1.1 is the pair shown 
in Fig.^which has one of the smallest impact parameter of our 
pairs. Provided the lensing galaxy is an elliptical, the magnifi- 
cation factor for this QSO is fi ~ 1.3. The system has an un- 
usual morphology and we cannot exclude the possibility that 
a more careful lensing analysis might give more interesting 
results. However, such detailed investigations are outside the 
scope of this paper. 

6. Summary and conclusions 

We have investigated the possibility to detect the reddening 
of QSOs by foreground galaxies and to measure the dust ex- 
tinction properties. Cross-correlating coordinates of the SDSS 
quasar catalogue with the New York University Value-Added 
Galaxy Catalog we find 164 QSO-galaxy pairs with impact pa- 
rameters within 30 kpc at the redshift of the foreground galaxy. 
With additional colour cuts we identified two low-z systems 



Redshift 
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Fig. 7. The best fit values of Ry ordered after galaxy redshift 
for the CCM and Fitzpatrick parametrisations. 

(z < 0.08) involving spiral galaxies where the QSO colours are 
particularly indicative of reddening in the intervening galaxy. 
Both QSOs have very large impact parameters in terms of the 
the optical radii of the foreground galaxies and are thus indica- 
tive of large dust column densities outside the optical disc of 
spiral galaxies. With this assumption, and using redshifted red- 
dened QSO templates we fit for Ry and E(B - V) using the 
measured colours and redshifts of the quasar and galaxy and 
find that Galactic type dust provides a good fit to the data. 
Such an interpretation needs to be confirmed by e.g. spectro- 
scopic means. We emphasize that we cannot rule out other ex- 
planations for the red colour of these two quasars, e.g. intrinsic 
redder colour, extinction in the host galaxy or in an unresolved 
foreground galaxy. 

We also explored the extinction properties in five QSOs in 
the literature, where the foreground galaxies have redshifts in 
the range z = 0.15 - 1.5. A wide range of preferred values of 
Rv was found. This is an indication that the dust properties in 
other galaxies may be different from the Milky Way. 

We note with interest that the dust in galaxies with high red- 
shifts in our study were better fitted by lower values of Ry than 
the average Milky Way value (See Fig. [7}. If this effect is con- 
firmed with a larger sample, and selection effects can be ruled 
out, it would be indicative of an evolution of the dust properties 
with redshift, which would have an impact on e.g. the use of ex- 
tinction corrected supernovae as distance indicators. It would 
not be too surprising if such an evolution in the reddening law 
exists since e.g. metallicity and elemental abundance ratios are 
known to change with redshift and these are expected to affect 
the dust properties. However, it is possible that the differences 
between our two sets are mainly due to selection effects since 
the few systems studied have been chosen with different meth- 
ods: at high-z through detection of absorption features and at 
low-z by selecting reddened QSOs at small impact parameters 
with respect to resolved foreground galaxies. 
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